Cosmological Constraints on f(R) Acceleration Models 
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Models which accelerate the expansion of the universe through the addition of a function of 
the Ricci scalar f(R) leave a characteristic signature in the large-scale structure of the universe 
at the Compton wavelength scale of the extra scalar degree of freedom. We search for such a 
signature in current cosmological data sets: the WMAP cosmic microwave background (CMB) power 
spectrum, SNLS supernovae distance measures, the SDSS luminous red galaxy power spectrum, 
and galaxy-CMB angular correlations. Due to theoretical uncertainties in the nonlinear evolution of 
f(R) models, the galaxy power spectrum conservatively yields only weak constraints on the models 
despite the strong predicted signature in the linear matter power spectrum. Currently the tightest 
constraints involve the modification to the integrated Sachs- Wolfe effect from growth of gravitational 
potentials during the acceleration epoch. This effect is manifest for large Compton wavelengths in 
enhanced low multipole power in the CMB and anti-correlation between the CMB and tracers of 
the potential. They place a bound on the Compton wavelength of the field be less than of order the 
Hubble scale. 

PACS numbers: 



I. INTRODUCTION 

Cosmic acceleration can arise either from an unknown 
component of dark energy with negative pressure or a 
modification to gravity that only appears at cosmological 
scales and densities. Additional terms in the Einstein- 
Hilbert action that are non- linear functions f(R) of the 
Ricci scalar R have long been known to cause acceleration 
P, H, d, 13 and have been the subject of much recent 




IS EE IE EE 

|2j,l25|,l26|,l27|,|28|,|2l,[33. 

The main challenge for f(R) models as a complete de- 
scription of gravity lies with the extremely tight con- 
straints on such modifications placed by solar system 
and local tests of general relativity. Chiba [3l[ showed 
that the fundamental problem is that f(R) models gener- 
ically introduce a light scalar degree of freedom with 
a long Compton wavelength at cosmological densities 
[32l [33 | . This problem can be mitigated by the so-called 
chameleon mechanism [34], [35[ , where the local Compton 
wavelength can decrease in high density environments 
[H, [3?J • The cosmological Compton wavelength can then 
be a scale of cosmological interest. Nevertheless, it must 
typically be less than a few tens of Mpc if the Galactic 
gravitational potential is not substantially deeper than 
implied by local rotation curve measurements [381 ] . 

In this Paper, we take the perspective that f(R) mod- 
els are effective theories that are valid at cosmological 
densities and scales and arc not necessarily predictive at 



the high densities and small scales of local tests. Hence 
we explore models with Compton wavelengths out to the 
horizon size and seek to constrain them from cosmologi- 
cal observables alone. At the very least, this exploration 
yields a cosmological test of general relativity that is in- 
dependent of local constraints. 

Although f(R) models can change the expansion his- 
tory during the acceleration epoch, there in principle al- 
ways exists a dark energy model that provides the same 
history. The unique and strongest signatures of f(R) 
modifications are on cosmological structure formation 
[25| and can impact observables down to Compton wave- 
lengths of a few Mpc [38[ . We consider constraints aris- 
ing from the CMB angular power spectrum measured by 
WMAP [3^|, the linear matter power spectrum inferred 
from the Sloan Digital Sky Survey (SDSS) luminous red 
galaxy (LRG) sample [4fJ, the distance measures by the 
Supernovae Legacy Survey (SNLS) [4l|, and the cross- 
correlation between the CMB and large scale structure as 
measured by WMAP and a range of galaxy and quasar 
surveys [H EE 0, HE HI . 

The outline of the paper is as follows. In ijTTl we review 
the f(R) model and discuss our calculation and analysis 
methods. In ijllll we present the results of joint cosmolog- 
ical constraints on these models. We discuss these results 
in JEVj 



II. f{R) METHODOLOGY 
A. Model 
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In f(R) models of gravity, the Einstein-Hilbert action 
is supplemented by a term that is non-linear in the Ricci 
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scalar R 



S 



R + f(R) 
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where C m is the matter Lagrangian. The modified Ein- 
stein and Friedmann equations result from varying the 
action with respect to the metric. Given the freedom to 
choose a functional form for f(R\ any desired expan- 
sion history can be replicated @, [2l Ha S3, H 51 ■ In 
particular, one can choose the ACDM expansion history 
which is known to satisfy distance constraints from high 
redshift supernovae, baryon acoustic oscillations and the 
CMB. 

Even given the degeneracy with dark energy in the 
expansion history, f(R) models have distinguishable ef- 
fects on the formation of structure. The promotion of 
the Ricci scalar to a dynamical degree of freedom modi- 
fies the force law between particles. This modified force is 
mediated by a new scalar degree of freedom Jr = df /dR, 
which has a squared Compton wavelength proportional 
to fun = d 2 f /dR 2 . Below the Compton wavelength scale 
gravity becomes a scalar-tensor theory, leading to an en- 
hancement in the growth of cosmological density pertur- 
bations and a corresponding suppression in the decay of 
gravitational potentials. 

For cosmological tests, it is convenient to express the 
squared Compton wavelength in the background in units 
of the Hubble length squared [H, HI 



B = 
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where a is the scale factor and the Ricci scalar is eval- 
uated at the background density. We will specialize our 
consideration to f(R) models that exactly reproduce the 
ACDM expansion history to test whether the unique sig- 
natures of f(R) gravity are seen in current cosmological 
data sets. 

Under this assumption, the additional degree of free- 
dom in f(R) gravity is parameterized by the value of 
B today, B a = B(z — 0). Stability requires the mass 
squared of the scalar to be positive (i.e. a prior of B n > 0) 
[H, [H. Note that in the limit that B -> 0, the phe- 
nomenology of ACDM is recovered in structure forma- 
tion tests as well as expansion history tests. More gener- 
ally, the control parameter is the average Compton wave- 
length through the acceleration epoch when gravity is 
modified. 



B. Power Spectra Calculation 

The fundamental observables of our f(R) model are 
the same as in ACDM. These include the redshift- 
distance relation, the CMB angular power spectrum Ci, 
galaxy power spectra P g (k), and the angular correlation 
between galaxies and the CMB w(8). Cosmic shear, dark 
matter halo profiles and masses from weak lensing as 



well as the cluster abundance are also potential observ- 
ables but their utilization requires cosmological simula- 
tions that are beyond the scope of this work. 

We modified the Boltzmann code CAMB [HTj] to cal- 
culate these observables in f(R) gravity. In the CAMB 
code, the density perturbations are computed in the syn- 
chronous gauge. We evolve the usual Boltzmann code up 
to a ~ 0.01 when deviations introduced by f(R) are still 
small. At this epoch, we transform the matter pertur- 
bations from synchronous gauge to gauge invariant vari- 
ables, and feed them as initial conditions to the linear 
perturbation equations for the density fluctuation and 
CMB sources (see [25| Eqs. 28-35). Once the initial con- 
ditions are specified, the power spectrum observables are 
computed from a separate f(R) routine that bypasses 
the usual CAMB code without significantly increasing 
the computational time. 

As we shall discuss, the main modification made by 
our f(R) models on the CMB is a change in the evolution 
of gravitational potentials during the acceleration epoch. 
The CMB power spectrum is modified at low multipolcs 
through the so called "Integrated Sachs- Wolfe" (ISW) ef- 
fect, and the temperature field is correlated with tracers 
of gravitational potentials such as galaxies. This corre- 
lation is strong if the redshift of the galaxies is matched 
to the epoch at which the gravitational potentials evolve. 
Following [52| , we model the angular correlation between 
a set of galaxy surveys indexed by i and the CMB, as- 
suming a broad band selection function for the galaxies, 
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(3) 



where the Zi are chosen to match the median redshift in 
the surveys. While crude, the uncertainties introduced 
by this choice of selection function are smaller than the 
difference in model predictions that will be considered. 
The galaxy bias is assumed to be constant in each red- 
shift bin. We calculate the cross power spectrum between 
gala xy number density and CMB temperature following 
|25| (Eqs. 55-60) and transform it to the angular power 
correlation function Wi(9). 



C. Likelihood Analysis 

We use a Markov Chain Monte Carlo (MCMC) tech- 
nique [H, [H, [El, [H, [13, [El to evaluate the likeli- 
hood function of model parameters. The MCMC is used 
to simulate observations from the posterior distribution 
V{0\x), of a set of parameters 9 given event x, obtained 
via Bayes' Theorem, 



V(9\x) = 



V(x\9)V(9) 
fV(x\9)V(9)d0' 



(4) 



where V(x\6) is the likelihood of event x given the model 
parameters 9 and V(0) is the prior probability density. 
The MCMC generates random draws (i.e. simulations) 
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from the posterior distribution that are a "fair" sample 
of the likelihood surface. From this sample, we can es- 
timate all of the quantities of interest about the poste- 
rior distribution (mean, variance, confidence levels). A 
properly derived and implemented MCMC draws from 
the joint posterior density V{0\x) once it has converged 
to the stationary distribution. We use 16 chains and a 
conservative Gelman-Rubin convergence criterion [59| to 
determine when the chains have converged to the station- 
ary distribution. 

For our application, 6 denotes a set of cosmological 
parameters. We then use a modified version of the Cos- 
moMC code (5|| to determine constraints placed on this 
parameter space by WMAP, SDSS LRG galaxy power 
spectrum, and supernova distance measures. For the 
LRG data, we only use the first 14 fc-bins in our fits 
reflecting a conservative cut for linearity k < 0.1/i/Mpc. 
For the supernovae constraint, we take the Supernovae 
Legacy Survey data set and its analysis remains unaf- 
fected by our f(R) modification. Likewise, its impact is 
mainly to help determine expansion history parameters. 
In our analysis, we take the parameter set {tub = f^/i 2 , 
u} m = fl m h 2 , e A ,Ml0 lt) A s ) 1 n s ,T 1 B Q ,b,Q ni }, where 6 A 
is the angular size of the acoustic horizon, and A s is 
the power in the primordial curvature perturbation at 
k = 0.05 Mpc" 1 . The universe is assumed to be to be 
spatially flat. Recall also that we fix the expansion his- 
tory to be the same as a ACDM model. Hence 9a is also 
a proxy for the Hubble constant Ho = lOO/i km/s/Mpc 

or n m . 

The last two parameters in the set are the linear galaxy 
bias b and non-linearity parameter Q n \ required for the 
interpretation of the galaxy power spectrum P g {k) of 
SDSS LRGs [H (see also flTTBjl . defined as 



P g {k) = P L (k) b- 



1 + Qnifc 2 
1 + 1.4& 



(5) 



The linear bias factor b is defined with respect to the 
linear matter power spectrum Pi,(k) at z = 0.35, the ef- 
fective redshift of the LRG galaxies. Since f(R) models 
predict a scale dependent growth rate that changes the 
shape of the power spectrum as a function of redshift, 
we do not define the bias as relative to the power spec- 
trum at z — (c/. 40]). The second factor on the right 
hand side accounts for the non-linear evolution of the 
matter power spectrum shape and the scale-dependent 
bias of the galaxies relative to dark matter. We assume 
that this ansatz continues to be a valid approximation 
for f(R) theories. We do not marginalize analytically 
over the b and Q n i parameters as is done normally in the 
CosmoMC code; instead the parameters are marginal- 
ized numerically as independent nuisance parameters in 
the MCMC analysis. 

For all of our parameters except for Bq, we employ 
fiat linear priors in the stated parameters. For Bo, we 
supplement the flat prior with the stability condition that 
B > 0. 

Finally, these constraints are projected onto the space 



of the galaxy-ISW angular correlation function for com- 
parison with the data. We do not attempt here to include 
the correlation function data in the likelihood. That 
would require a joint re-analysis of the various data sets 
to capture the covariance between the measurements. It 
would also require an assessment of systematic errors 
in the correlation due to uncertainties in the selection 
functions and other effects. Instead we compare predic- 
tions with the individual measurements and their errors 
as quoted in the literature. 
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FIG. 1: Time evolution of the effective gravitational potential 
$_ = (*- *)/2 for k = 10" 3 /Mpc near the peak contribution 
to the ISW effect at low multipoles. As the Compton wavelength 
parameter Bq increases, the decay of the gravitational potential 
in ACDM decreases and eventually turns to growth. This reversal 
changes the sign of the ISW effect in overdense regions traced by 
galaxies. 



III. COSMOLOGICAL CONSTRAINTS 
A. CMB 

In f(R) models that follow the ACDM expansion his- 
tory considered here, none of the CMB phenomenology at 
recombination is affected by the modification to gravity. 
The Compton wavelength parameter B is driven rapidly 
to zero at high density and curvature. Hence all of the 
successes of ACDM in explaining the acoustic peaks car- 
ries over to these f(R) models. 

The impact of f(R) gravity on the CMB comes ex- 
clusively through the so-called Integrated Sachs- Wolfe 
(ISW) effect at the lowest multipole moments. The ISW 
effect arises from an imbalance between the blueshift a 
CMB photon suffers while falling into a gravitational 
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FIG. 2: CMB angular power spectrum G& for f(R) models with 
the Compton wavelength parameter Bo = (ACDM), 0.5, 1.5, 3.0, 
5.0. As Bo increases, the ISW contributions to the low multipolcs 
decrease, change sign, and then increase. WMAP3 data with noise 
error bars are overplotted and rule out Bq > 4.3 (95% CL). 



FIG. 3: Posterior probability distribution for the Compton wave- 
length parameter Bo inferred from the joint likelihood analysis of 
WMAP CMB angular power spectrum, SDSS LRG galaxy power 
spectrum, and SNLS supernovae data. The upper limit on Bq from 
the joint constraint is driven by the CMB data, specifically the ISW 
effect. 



potential well and the redshift while climbing out if 
the gravitational potential evolves during transit. With 
a cosmological constant, gravitational potentials decay 
during the acceleration epoch. For f(R) models, the en- 
hancement of the growth rate below the Compton scale 
can change the decay into growth (see Fig. [TJ. This re- 
versal changes the sign of the ISW effect. CMB photons 
then become colder along directions associated with over- 
dense regions. 

Fig. [5] illustrates the impact of this effect on the 
CMB power spectrum. As the Compton wavelength ap- 
proaches the Gpc (k ~ 10~ 3 Mpc^ 1 ) scales associated 
with the low multipoles of the ISW effect, the reduction 
in the decay of the potential also suppresses the ISW ef- 
fect. Near Bq w 1.5 the ISW effect is almost entirely 
absent at the quadrupole. Reduction in the amplitude of 
the quadrupole is in fact weakly favored by the data but 
the large cosmic variance of the low multipoles prevents 
this from being a significant improvement. Moreover, the 
elimination of the ISW effect at higher multipoles where 
the observed power is higher counteracts this improve- 
ment. For Bq > 1.5 the Compton wavelength exceeds 
the scales of interest and potential decay turns to po- 
tential growth. By Bo ~ 3 the ISW effect has an equal 
amplitude to that of ACDM. For Bo > 5, it is too large 
to accommodate the WMAP data. These features drive 
the overall joint constraint on Bq shown in Fig. [3J The 
WMAP data also serve to fix the parameters that control 
the high redshift universe. In particular, it constrains the 



initial amplitude of power on scales that are observed in 
galaxy surveys such as the SDSS LRG survey. 

B. Galaxy Power Spectrum 

The time-dependent Compton wavelength of f(R) 
models induces a more dramatic effect in the matter 
power spectrum during the acceleration epoch. In fact, 
average Compton wavelengths down to a few Mpc are 
potentially observable in the linear power spectrum. Un- 
der this scale, the enhanced growth rate leads to excess 
power relative to the same initial power spectrum deter- 
mined by the WMAP data. Unfortunately the overall 
change in power is degenerate with the unknown galaxy 
bias. However if the Compton wavelength appears be- 
tween the few to 100 Mpc scale the distortion of the shape 
of the power spectrum is potentially distinguishable in 
current surveys. In our parameterization this occurs for 
B ~ 10~ 5 - 10~ 2 . 

The SDSS LRG data set in fact favor enhanced power 
over the linear ACDM power spectrum. In Fig. [4] 
we compare a linear ACDM power spectrum (w& = 
0.025, w m = 0.128,^0 = 73,10*% = 21.2, n s = 
0.95, t = 0.09, B = 0,6 = 2.2 and a linear f(R) power 
spectrum with the same parameters and Bo = 4 x 10~ 4 . 
The f(R) model is in fact a better fit to the shape of the 
power spectrum, with Axis = 2AlnL « 11.5 for this 
specific choice of bias. 
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FIG. 4: SDSS LRG galaxy power spectrum data compared with a 
linear ACDM power spectrum (see text), a linear f(R) model with 
Bo = 4 X 10 -4 and the best fit nonlinear ACDM power spectrum. 
While the data are better fit by f(R) with linear power spectra, 
nonlinearities make ACDM a better fit to the data and open up a 
degeneracy between f(R) modifications parameterized by Bo and 
the non-linearity parameter Q n i. 



However in the ACDM model, we expect LRG galaxy 
clustering to be nonlinear in exactly the region where 
these changes of shape occur. In fact a non-linearity 
parameter of Q n \ — 30, which represents a reasonable 
amount of non-linearity, produces an excellent fit to the 
data (see Fig. |4j. Hence, f(R) enhancements of small 
scale power are degenerate with non-linear effects in 
ACDM and open up a corresponding degeneracy between 
Q n \ and Bq (see Fig. [5]). The non-linear modification also 
introduces a small suppression of power at intermediate 
k for any value of Q n \ which also marginally improves the 
fit. 

Furthermore if Bq > 10~ 2 the Compton scale exceeds 
the largest scales in the survey at k ~ 10~ 2 /i/Mpc. The 
enhancement of small scale power, while pronounced, be- 
comes degenerate with the galaxy bias b. With the pri- 
mordial amplitude fixed by WMAP, the bias must de- 
crease as Bq increases, leading to an anti-correlation be- 
tween the two parameters (see Fig. [6]). 

Even order unity Bq is allowed after marginalization. 
The LRG data do weakly disfavor even larger Bq but the 
overall constraint is dominated by the CMB data (see 
Fig. [3] and compare Bq = 3 relative to Bq = 0, where the 
CMB spectra are nearly identical). 

It is the marginalization over the non-linear parame- 
ter Q n \ and bias which substantially degrades the abil- 
ity of the LRG data to constrain f(R) models. This 
is a theoretical and not an observational limitation. As 
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FIG. 5: The 2D joint 68% CL (dark) and 95% CL (light) 
constraints on the Compton wavelength parameter Bq and non- 
linearity parameter Q n \ inferred from the joint likelihood analysis 
of WMAP CMB angular power spectrum, SDSS LRG galaxy power 
spectrum, and SNLS supernovae data. The upper panel focuses in 
on the range Bq < 0.005 in order to highlight the anti-correlation 
between these two parameters discussed in the text, while the lower 
panel shows the full joint constraint. 



we have seen, even the small value of Bq ~ 10~ 4 sub- 
stantially impacts the current data. With cosmological 
simulations that address the non-linear evolution of the 
matter power spectrum in f(R) and the association of 
LRG's with dark matter haloes in the simulation, this 
uncertainty can be lifted leading to substantially tighter 
bounds on the Compton wavelength. 



C. Galaxy-ISW Correlation 

The angular correlation between the CMB tempera- 
ture field and galaxy number density field induced by 
the ISW effect places an interesting constraint on f(R) 
models with averaged Compton wavelengths in the > 100 
Mpc regime. As discussed in IIII A[ the ISW effect re- 
verses sign for wavelengths smaller than the average 
Compton wavelength during acceleration. This rever- 
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FIG. 6: The 2D joint 68% CL (dark) and 95% CL (light) con- 
straints on the Compton wavelength parameter Bo and the linear 
galaxy bias b inferred from the joint likelihood analysis of WMAP 
CMB angular power spectrum, SDSS LRG galaxy power spectrum, 
and SNLS supernovae data. The upper panel focuses in on the 
range Bo < 0.005 in order to highlight the anti-correlation be- 
tween these two parameters discussed in the text, while the lower 
panel shows the full joint constraint. 



sal of sign causes galaxies to be anti-correlated with the 
CMB [2a| . For sufficiently large Bq > 1, the Compton 
wavelength subtends a scale larger than the 1 — 10° over 
which the correlation has been measured. In this regime, 
f(R) models predict that galaxies are anti-correlated 
with the CMB, in conflict with the data measure d by 
2MASS, APM, SDSS, NVSS and QSO [I| SB El, [If. 

Figure [7] shows the angular correlation data at 6° from 
a compilation in [&oj ] (updated by Gaztanaga, private 
communication). We compare these data with predic- 
tions from models in the chain spanning the 68% and 
95% allowed regions given by the WMAP CMB, SDSS 
LRG power spectrum and supernovae data sets. The 
lack of anti-correlation in any given data point places an 
upper bound of Bq < 1 at the significance level of the 
detection. This is approximately a factor of 4 improve- 
ment in Bq over the other constraints or a factor of 2 in 
the Compton wavelength. These individual constraints 



redshift z 

FIG. 7: The galaxy-ISW angular correlation function at 6 degrees 
is shown as a function of redshift as measured by a variety of surveys 
(including quasars), with the bias divided out. The constraints on 
f(R) models given by CMB, SDSS LRG galaxies and supernovae 
are projected onto this observable, with 68% CL (dark) and 95% 
CL (light) shaded regions shown. The ACDM limit of Bq — > is 
shown by a solid line, and the dashed lines (from top to bottom) 
correspond to So = {0.5, 1.5, 3.0, 5.0}. 



can be improved somewhat by a full joint analysis of the 
correlation measurements but that lies beyond the scope 
of this work. 



IV. DISCUSSION 

We have analyzed current cosmological constraints on 
f(R) acceleration models from the CMB, SDSS galaxy 
power spectrum, galaxy-ISW correlations and super- 
novae. By choosing f(R) models with a ACDM expan- 
sion history, we have explored whether the unique sig- 
natures of the f(R) modification to gravity are seen in 
current data. 

Despite the relatively large impact of f(R) models on 
the matter power spectrum, the strongest current con- 
straints involve the modification these models induce on 
the ISW effect in the CMB. The growth of density per- 
turbations is enhanced under the Compton wavelength 
of the field fn = df / dR by scalar-tensor modifications to 
the gravitational force law. This enhancement turns the 
decay of gravitational potentials during the acceleration 
epoch in ACDM into growth. The joint constraint on 
the Compton wavelength parameter Bq < 4.3 (95% CL), 
from the WMAP CMB power spectrum, SDSS galaxy 
power spectrum and supernovae, is driven by the CMB. 
This constraint allows Compton wavelengths as large as 
the current horizon and is a weak bound on f(R) models. 

A stronger bound of Bq < 1 can be inferred from the 
positive correlation between the CMB and a range of 
galaxy surveys from z ~ 0.1 — 1.5. In ACDM this positive 
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correlation is induced by the ISW effect from a decaying 
potential. In f(R) models growing potentials convert 
the positive correlation into anti-correlation in violation 
of the observations. 

The weak impact of the SDSS LRG galaxy power spec- 
trum on our joint constraints is due to a theoretical 
rather than observational limitation. Intriguingly the lin- 
ear f(R) power spectrum is a marginally better fit to 
the data than a linear ACDM one. However the nonlin- 
ear corrections expected in ACDM bring the model back 
in agreement with the data. Y-body simulations have 
yet to determine the non-linear matter power spectrum 
and dark matter halo (and hence galaxy) correlations in 
f(R) models. The likely impact of non- linear evolution 
is degenerate with the enhancement of small scale linear 
power seen in f(R) models. When non-linearity of the 
form expected for ACDM is marginalized for f(R) mod- 
els, the impact of the galaxy power spectrum on joint con- 
straints becomes very weak, and Compton wavelengths 
out to the horizon length are allowed. 

In the future, the window between the horizon length 



and a few tens of Mpc can be tested by larger photomet- 
ric surveys that probe the shape of the power spectrum 
across the whole range of scales. In principle, Compton 
wavelengths down to a few Mpc can be tested by galaxy 
power spectra and cosmic shear from comparison of data 
with simulations of f(R) models. 
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